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Introduction: In this work, we will compare ob- 
served temperatures of the surface of Vesta using data 
acquired by the Dawn [1] Visible and Infrared Map- 
ping Spectrometer (VIR-MS) [2] during the approach 
phase to model results from the KRC thermal model. 

Thermal Inertia: Thermal inertia (I) is a measure 
of how quickly a material responds to changes in tem- 
perature and is defined as 1=^1 kpc, where k is the ther- 
mal conductivity, p is the bulk density, and c is the 
specific heat capacity. SI units for / are J m" 2 s" 1/2 K" 1 
which henceforth will be refer to as “TIU” (thermal 
inertia unit) as proposed by Putzig [3]. 

High thermal inertia materials, such as bedrock, re- 
sist changes in temperature while temperatures of low 
thermal inertia material, such as dust, respond quickly 
to changes in solar insolation. The surface of Vesta is 
expected to have low to medium thermal inertia val- 
ues, with the most commonly used value being ex- 
tremely low at 15 TIU [4]. 

There are several parameters which affect observed 
temperatures in addition to thermal inertia: bond al- 
bedo, slope, and surface roughness. In addition to 
these parameters, real surfaces are rarely uniform 
monoliths that can be described by a single thermal 
inertia value. Real surfaces are often vertically layered 
or are mixtures of dust and rock. For Vesta’s surface, 
with temperature extremes ranging from 50 K to 275 K 
and no atmosphere, even a uniform monolithic surface 
may have non-uniform thermal inertia due to tempera- 
ture dependent thermal conductivity. 

KRC Thermal Model: We employ a multi- 

layered thermal-diffusion model called ‘KRC’ [5], 
which has been used extensively in the study of Mar- 
tian thermal physical properties [e.g. 5-11]. This ther- 
mal model is easily modified for use with Vesta by 
replacing the martian ephemeris input with the Vesta 
ephemeris and turning the atmosphere off. This model 
calculates surface temperatures throughout an entire 


Vesta year for a specific set of slope, azimuth, latitude, 
elevation and albedo for up to ten different thermal 
inertia values. To calculate surface temperatures, KRC 
uses the Delta-Eddington approximation for radiative 
flux (Kieffer et al., 1977) to solve the subsurface ther- 
mal diffusion equation using finite-difference methods. 
The upper boundary condition is solar insolation (cal- 
culated at each step according to orbital position, or- 
bital inclination and time of day). KRC models are run 
with varying values for albedo, thermal inertia and 
layer depth to create temperature indices which are 
closely matched, via bilinear interpolation, to the dates 
and times observed by VIR. After this interpolation, 
our working index contains modeled temperatures for 
all times of observation at a logarithmic range of ten 
thermal inertias between 1 5 and 400 TIU. Higher val- 
ues of TIU can be included if necessary. The downhill- 
simplex method of function minimization (Nelder and 
Mead, 1986) is then used to compare the observed and 
modeled data to determine best-fit values for thermal 
inertia (and layer depth when applicable). 

Vesta Approach Data: The Dawn spacecraft [1] 
acquired approach observations on June 30, 2011. 
Figure 1 shows a false-color view of Vesta as seen by 
VIR in the near-infrared. These data provide full disk 
view of Vesta, providing a range of local times (Fig 2). 



Figure 1: VIR IR False color image. 


Vesta Temperature Data: Temperatures were 

calculated using a Bayesian approach to nonlinear in- 
version as described by Tosi et al. [13, 14]. In order to 
compare observed temperatures of Vesta to model cal- 
culations, several geometric and photometric parame- 
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ters must be known or estimated. These include local 
mean solar time, latitude, local slope, bond bolometric 
albedo, and the effective emissivity at 5 pm. Local 
time, latitude, and local slope are calculated using the 
USGS ISIS [15, 16] software system. The bolometric 
bond albedo is estimated from the observed surface 
reflectivity. The effective emissivity is estimated by 
comparing the total radiance at 5 pm to the solar flux 
and the thermal emission. If one assumes that Kirch- 
hoff s Law applies and that the surface reflection can 
be approximated as Lambertian, then the observed 
radiance (R G b s ) is a function of emissivity: 


R obs = (1-A^iio so i+£)fix (T)=S X (ju 0 sd +Px (T)) 
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Figure 2: Comparison of observed temperatures to 
KRC model results assuming a smooth flat surface. 
The plus signs are VIR temperatures and the 
dashed lines are the KRC models. The colors indi- 
cate latitude where white is 80°S, red is the latitude 
range 75°S to 20°S, green is the latitude range 20°N 
to 20°S and blue is the latitude range 20°N to 60°N. 

Roughness Effects: Because KRC is a 1-D model, 
the effects of roughness are approximated by combin- 
ing a suite of KRC model results with a range of 
slopes and slope azimuths. While this approach cap- 
tures many of the effects of surface roughness, it does 
not correct for the secondary effects of radiative cou- 
pling between facing surfaces of varying temperatures. 

Results: Figure 2 compares the observed tempera- 
tures to KRC model results assuming a smooth flat 
surface with a fixed bond albedo of 0.2 and a surface 
emissivity of 0.6. The model curves shown are for a 
thermal inertia of 30 TIU. The need to incorporate 
surface roughness is clearly illustrated. The tempera- 
ture profiles as a function of latitude and local time of 
day are generally the correct shape, but effects from 
beaming and surface roughness are needed to properly 
model the warm temperatures observed. 

Future Work: The incorporation of local slopes, 
emissivity, and albedo is the next step in this project. 

Conclusion: Thermal inertia analysis of the sur- 
face of Vesta has great potential to identify a variety of 
processes ranging from variations in grain size to ef- 
fects of cementation. 



